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Abstract We examined the applicability of end-notched 
flexure (ENF) tests for measuring the mode II fracture 
toughness of wood. Western hemlock (Tsuga heterophylla 
Sarg.) was used for the specimens. The fracture toughness 
at the beginning of crack propagation G !Ic and that during 
crack propagation G IIR were calculated from the load¬ 
loading point compliance and load-crack shear displace¬ 
ment (CSD) relations. The obtained results were compared 
with each other, and the validity of measurement methods 
were examined. The results are summarized as follows: (1) 
The value of G nc increased with the increase in initial crack 
length. When measuring G llc by ENF tests, we should be 
aware of the dependence of G„ c on the initial crack length. 
(2) The value of G IIR initially increased with the crack 
length, and it reached a constant value. (3) Measurement 
of the CSD is recommended when obtaining G IIR because 
the crack length, which has a great influence on the G UR 
calculation, is implicitly included in the CSD. (4) We found 
that the crack length during its propagation should be 
evaluated by the final crack length. 
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Introduction 


In linear-fracture mechanics, an in-plane shear mode 
fracture is defined as a “mode II” fracture; and the mode II 
fracture behaviors of wood had been variously analyzed. 17 
Recently, a draft standard for measuring the mode II 
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fracture toughness was proposed by the International 
Union of Testing and Research Laboratories for Materials 
and Structures (RILEM). 7 Nevertheless, the specimen 
shape in the draft is complicated, and we are afraid that the 
proposed method cannot become popular. 

Among the various mode II fracture testing methods, 
the end notched flexure (ENF) test is frequently used for 
fiber-reinforced composites, and testing techniques have 
been developed during the past two decades. 8 ' 13 We thought 
that these experimental techniques might be applied to 
measurement of the mode II fracture toughness of wood. 
Here, we examined the validity of the ENF testing method 
for measuring the mode II fracture toughness of wood. 


Theoiy 

Energy release rate derived by ENF test 

Russell and Street derived an equation that gives the mode 
II energy release rate by ENF test. 8 As shown in Fig. 1, an 
ENF specimen with the crack length of a is loaded by P at 
the center of the span of 2 L. In this case, the equation of 
bending is given by the elementary beam theory as follows. 


E X I d 2 y 1 /. % 

—-4- = — Px 0 <x<a) 

8 dx 2 4 v ’ 

■ E X I^~ = ——Px la^x^L) 

x dx 2 2 v ’ 

E x I^r = -Px-PL (L<jc<2L) 

dx 2 2 v ’ 


( 1 ) 


where E x is Young’s modulus in the longitudinal axis, and I 
is the second moment of cross-sectional area in the crack- 
free region. Solving this equation, we can obtain the dis¬ 
placement at the loading point, v, as follows. 

Pills' + 3a 3 ) 
v = —-- - 
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The loading point compliance C is defined as follows: 


c= — = 2L3+3fl3 
~ P ~ 12 E X I 

The energy release rate, G n , is given as: 


(3) 


= 

2b da 


(4) 


where b is the width of the beam. By substituting Eq. (3) 
into Eq. (4), the energy release rate can be written as: 


, = 9CP 2 a 2 
r “ 2b(2L 3 + 3a 3 ) 


(5) 


Another equation for the energy release rate was given 
by Kageyama and colleagues. 10 The crack shear displace¬ 
ment (CSD) d, which is defined in Fig. 1, is obtained by 
integrating the strain along the longitudinal axis from the 
lower end corner of the upper beam to the crack tip; it is 
given as follows. 
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where h is the upper beam height of cracked region, /' is 
the second moment of cross-sectional area in the cracked 
region (which equals 1/8), and M is the bending moment. 
The compliance A in this case is written as follows: 


6 _ ha 2 
T ~ 8 EJ' 


(7) 


From Eqs. (3), (5), and (7), the energy release rate can be 
formulated without the crack length as: 


P 2 dC _ P 2 3 a 2 = 3 Pd = 3 P 2 l 
2b da ~ 2b 4E X I 8 ' bh 8 bh 


( 8 ) 


Stability of crack propagation 

The stability of crack propagation is evaluated by the sign 
of dGJda. Under the constant loading point displacement 


condition, dG u /da is derived from Eqs. (4) and (6) as 

follows. 10 


dG u _ 3v 2 a f ^ _ 9 a 3 'j 
Ha ~ \EJbC 2 { 2L 3 + 3a 3 j 


When the value of dG n lda is negative, the crack propaga¬ 
tion is stable. Then the fracture toughness during the crack 
propagation can be obtained when the crack length satisfies 
the following relation. 



= 0.7L 


( 10 ) 


Otherwise, the crack grows unstably, and only the fracture 
toughness at the beginning of crack propagation can be 
obtained. 


Experiment 

Materials 

Western hemlock (Tsuga heterophylla Sarg.) whose density 
of 0.48 g/cm 3 was used for the specimens. Specimens 
were conditioned at 20°C and 65% relative humidity (RH) 
before and during the tests. 


ENF tests 

A beam specimen was cut with the dimensions of 20mm 
(radial direction) X 20 mm (tangential direction) X 500mm 
(longitudinal direction). The crack was produced in the 
longitudinal-tangential plane, which is the so-called TL 
system. 1 The crack, with the thickness of 1 mm, was first cut 
by a band saw and then was extended by a razor blade. The 
initial crack length a 0 varied from 69 to 207 mm at intervals 
of 23 mm. This range corresponded to the value of a<JL from 
0.3 to 0.9 at intervals of 0.1. Two sheets of Teflon film with 
0.5 mm thickness were inserted between the crack surfaces 
to reduce the friction between the upper and lower cantile¬ 
ver beams. This specimen was supported with the span of 
460 mm, and a load was applied at the center of the spans at 
the crosshead speed of 2mm/min. 

The loading point displacement v was measured by a dial 
gauge set below the loading nose, whereas the crack shear 
displacement d was measured by the CSD gauge set on the 
upper cantilever of the specimen as in Fig. 2. As shown 
in Fig. 3, the “critical load” P c was determined as the load 
of the intersection point between the two straight-line 
segments through the pre- and postlinear portions of the 
load-loading point displacement or load-CSD curves. 

The fracture toughness at the beginning of crack pro¬ 
pagation C? tic was calculated for every specimen from Eqs. 
(5) and (8). For the specimen whose aJL was 0.7 or 0.8, the 
fracture toughness during the crack propagation G IlR was 
also obtained because the crack propagation was stable and 
the crack path was long enough in this initial crack length 
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Fig. 2. Crack shear displacement (CSD) gauge 



Loading point displacement v or CSD <5 


Fig. 3. Determination of the critical load P c 


range. Prior to the tests, straight lines were drawn in the 
crack-free region of an LT plane at intervals of 5 mm per¬ 
pendicular to the crack. A load was applied until a marker 
line split. This process was performed under a magnifying 
glass. The specimen was then unloaded until a certain load 
was attained, and the crack length was measured again 
because the crack opened a little during the unloading. 
The specimen was then reloaded until another 5-mm crack 
appeared. This load-unload process was repeated, and 
the value of G IIR corresponding to each crack length was 
obtained. 

In ENF test, the crack propagates without a clear open¬ 
ing, and it is difficult to measure the propagating crack 
precisely. Therefore, we tried to estimate the crack length 
by the following five equations. From Eq. (5), the propagat¬ 
ing crack length a p is derived as: 8 


fl p = 


Co 3 2 

fc p ) 



- 1 

L 3 

C 0 3 

lCo J 



( 11 ) 


where a 0 and C„ are the initial crack length and the initial 
loading point compliance, respectively; and C p is the loading 
point compliance during crack propagation. Similarly, a p is 
derived by the final crack length a f and loading point com¬ 
pliance C f as follows. 



C„ , 2 

(C p ) 


a n = 

-S-4+- 

-E- ~ 1 

L 3 

p 

Cf 3 

[C, J 



From Eq. (7), a p is derived as: 



( 12 ) 


(13) 


where 1 0 and X p are the load-CSD compliances at the initial 
and during the crack propagation, respectively. Similarly, 
a p is derived by the final crack length a t and the CSD 
compliance X ( as: 


•’-“ft (14) 

Another equation, proposed by Kageyama and colleagues, 
is shown as follows. 12 
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(16) 


and 



(17) 


where E y is Young’s modulus in the direction of thickness, 
and G xy is the shear modulus. The values of E x , E , and G xy 
were determined by compression tests (see below). 


Compression tests 

To determine the values of y3 0 and />, in Eq. (15), we mea¬ 
sured Young’s moduli E x and E y and the shear modulus G xy 
by compression tests. Short-column specimens whose 
dimensions were 40 X 20 X 20 mm were prepared. When 
measuring E x and E y , the long axis of the specimen coin¬ 
cided with the longitudinal and tangential directions of 
wood, respectively; whereas when measuring G xy , the long 
axis coincided with the direction inclined at 45 degrees 
with respect to the grain. Strain gauges were bonded at the 
centers of the LT planes; and from the stress-strain relation, 
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Young’s moduli were obtained. The shear modulus were 
determined from the following equation. 


s~i _ ^45 

'-'xy ~ 


2(1 + V 45 ) 


(18) 


where E 45 and v 45 are Young’s modulus and Poisson’s ratio 
in the direction inclined at 45 degrees with respect to the 
grain, respectively. 


Results and discussion 

Fracture toughness at the beginning of crack propagation 

Figure 4 represents the fracture toughness at the beginning 
of crack propagation G llc corresponding to the initial crack 
length a 0 . The value of G llc increased with the increase of the 
initial crack length a 0 whether the compliance was deter¬ 
mined from the loading point displacement or CSD. Hence, 
we believe that the value of G llc should be reported with the 
initial crack length as an important testing condition. 


propagation length Aa. The fracture toughness of brittle 
material tends to be constant even during the early stage of 
crack propagation. 11 In contrast, it is reported that in mate¬ 
rials with irregular fiber distribution a bridging between the 
crack surfaces occurs; the crack propagation is interrupted 
by this “fiber bridging,” and the G IIR increases with crack 
propagation (K. Kageyama, personal communication, 
1999). In our experiment, the value of G IIR reached a con¬ 
stant value when the crack grew to more than 200 mm. The 
representative value of G IIR , which is denoted as G IIR , was 
then defined as the average of the data obtained from the 
>200mm crack length region. Table 1 shows the values of 
G iir obtained from the load-loading point compliances and 
load-CSD compliance relations. We think that the value of 
G iir varies by the estimation method, and that the method 
for determining G IIR should be examined furthermore. 
As mentioned by Tanaka and colleagues, it may be advan¬ 
tageous to evaluate fracturing behavior from the whole 
relation between fracture toughness and crack propagation 
length. 11 


Fracture toughness during the crack propagation 

Figure 5 shows examples of the load-loading point displace¬ 
ment and load-CSD relations. For the ENF tests of brittle 
epoxy material, a maximum load can be observed in one 
load-unload procedure. 8 In our experiment, the load did not 
decrease during crack propagation. 

Figure 6 presents examples of fracture toughness during 
the crack propagations G nR corresponding to the crack 


Table 1. Rep rese ntative value of fracture toughness during crack 
propagation, G lIR 


aJL 

Measurement of compliance (kJ/m 2 ) 


Loading point compliance 

Crack shear displacement 

0.7 

0.61 ± 0.13 

0.68 ± 0.11 

0.8 

0.67 ± 0.27 

0.71 ± 0.11 


Results are the average ± SD 



Initial crack length/half span a^/L 



0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 
initial crack length/half span a q/L 


Fig. 4. Fracture toughness at the beginning of crack propagation G IIc 
corresponding to the initial crack length/half span. Circles and horizon¬ 
tal bars represent the mean and standard deviations, respectively. Frac¬ 


ture toughness was calculated from the loading point compliance (left) 
and the CSD compliance (right) 
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Fig. 5. Examples of the records of the load-unload process. Initial crack length/half span was 0.7 




Crack propagation length Aa (mm) Crack propagation length Aa (mm) 

Fig. 6. Changes in the fracture toughness G lm during crack propagation. Open circles and solid circles were obtained from the loading point 
compliances and CSD compliances, respectively; dashed line represents the overall crack length of 200mm 


When calculating fracture toughness from the load-CSD 
relation, the crack length is implicitly included in the CSD 
compliance X, and we can obtain the value of G IIR without 
measuring the crack length. Considering the difficulty of 
monitoring the propagating crack length, CSD measure¬ 
ment is recommended for obtaining the fracture toughness 
G nR . Additionally, the crack can propagate stably when the 
CSD is controlled to increase at a constant rate. 1011,13 In our 
experiment, however, this “CSD control method” could not 
be used because our testing machine was incapable of such 
measurement. In the future, the “CSD control method” 
should be examined for obtaining details of the mode II 
fracture behaviors of wood. 


Measurement of crack length 

To evaluate the crack length by the equation proposed by 
Kageyama and colleagues, Eq. (15), we used values of 
12.3, 0.80, and 1.18 GPa for E„ E y , and G xy , respectively. 
Figure 7 shows a comparison between the calculated and 
observed crack lengths. The crack length calculated by Eqs. 

(11) and (13), which contains the initial crack length and 
compliances, deviated from the observed crack length even 
during an early stage. In contrast, when using the final crack 
length and compliances, the crack length calculated by Eqs. 

(12) and (14) were closer to the observed one when the 
crack length was long enough. The latter tendency was also 
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160 170 180 190 200 210 220 


Measured crack length a (mm) 

Fig. 7. Comparison of the predicted and measured crack lengths. Ini¬ 
tial crack length/half span was 0.7 


found in the results obtained from Eq. (15), proposed by 
Kageyama and colleagues; and the coincidence was good 
in the crack length range where G IIR was a constant value. 
Considering the difficulty of measuring the propagating 
crack length, we recommend that it be evaluated by means 
of the final crack length. 


Conclusions 

We measured the mode II fracture toughness of wood 
by end-notched flexure tests, and following results were 
obtained. 

1. The fracture toughness at the beginning of crack 
propagation G llc increased with the increase of the initial 
crack length. Hence, we should consider the dependence of 
the fracture toughness on the initial crack length. 

2. The fracture toughness during the crack propagation 
G iir initially increased and was close to a constant 
value. 


3. Measurement of CSD is recommended when obtain¬ 
ing G I1R because the crack length, which has a serious 
influence on the calculation of G IIR , is implicitly included in 
the CSD. 

4. We recommend that the crack length during propaga¬ 
tion be evaluated based on the final crack length. 
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